Entropy of thermal quasiparticles in nuclei 
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Information on level density for nuclei with mass numbers A ~ 20-250 is deduced from discrete 
low-lying levels and neutron resonance data. The odd-mass nuclei exhibit in general 4 — 7 times the 
level density found for its neighboring even-even nuclei at the same excitation energy. This excess 
corresponds to an entropy of ^ 1.7 fes for the odd particle. The value is approximately constant 
for all mid-shell nuclei and for all ground state spins. For these nuclei it is argued that the entropy 
scales with the number of quasiparticles. A simple model based on the canonical ensemble theory 
accounts qualitatively for the observed properties. 
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I. INTRODUCTION 



Nuclear structure at low excitation energy depends on 
the residual long-range two-body interaction. One of the 
most exciting consequences of this interaction is the form- 
ing of J = nucleon pairs, the so-called Cooper pairs, 
where nucleons are moving in time reversed orbitals. It 
has been a great challenge to observe the breaking of 
these Cooper pairs and the general quenching of pair 
correlations as function of temperature T and angular 
frequency m. Rotational high spin states are essentially 
described by quasiparticles coupled to the rotating un- 
derlying core. Backbending phenomena and spin align- 
ments have been described according to their quasiparti- 
cle properties. 

The concept of thermal quasiparticles is based on the 
idea that thermal properties at low excitation energy are 
governed by a few quasiparticles coupled to the cold core 
of Cooper pairs. These quasiparticles are thermally scat- 
tered on available single particle states. Thus, they are 
not well-defined by one single (Nilsson) orbital with given 
spin and parity, but exhibit an average of the spectro- 
scopical properties of the orbitals in the vicinity of the 
Fermi surface. At higher temperatures (T > 0.5 MeV), 
the pair correlations are quenched and the core of Cooper 
pairs is no longer well defined jl] . 

The entropy is a fundamental quantity in the descrip- 
tion of many body systems. For hot nuclei the entropy 
has the same importance as the spin alignment has for 
rotational nuclei. It can be evaluated from mechanical 
variables like energy, particle number and volume. Thus, 
it can be deduced in the microcanonical ensemble, which 
is the appropriate ensemble for an isolated system like 
the nucleus. The microcanonical entropy is defined for 
all isolated systems, from few-body systems up to infi- 
nite systems. It has been shown that phase transitions 
in small systems can be studied using the microcanonical 
entropy without invoking the thermodynamical limit . 



Although less appropriate, also the canonical ensem- 
ble theory is widely applied to nuclei. One advantage of 
the canonical ensemble is that the thermodynamical ob- 
servables are smooth functions of temperature due to the 
Laplace transformation involved in the calculation of the 
canonical partition function. Therefore, we will model 
our findings in this work within the canonical ensemble 
theory. 

In a recent paper j^, we extracted the microcanon- 
ical single quasiparticle entropy in hot i6i,i62]-)y 
i7i,i72y]-, nuclei to be 1.70(15) in units of the Boltzmann 
constant ks- Furthermore, the single quasiparticle en- 
tropy was found to be approximately constant in the 1- 
6 MeV excitation region. 

In this paper we will investigate the global systematics 
of nuclear entropy as function of nuclear mass, excitation 
energy (or temperature), number of excited quasiparti- 
cles and ground state spin. In Sect. II we sketch a simple 
model for thermal quasiparticles in nuclei. The method 
of determining level density by counting levels and using 
neutron resonance spacings is described in Sect. III. In 
Sect. IV we display the systematics of entropies and in 
Sect. V the properties of thermal quasiparticles are dis- 
cussed. In Sect. VI we present an application involving 
thermal quasiparticles and, finally, concluding remarks 
are given in Sect. VII. 



II. SIMPLE MODEL FOR THERMAL 
QUASIPARTICLES 

In this section we will describe a simple model demon- 
strating the properties of thermal quasiparticles, which 
in turn will be compared to experimental data. The main 
task is to describe the partition function Z in the canon- 
ical ensemble for a given temperature T and number of 
thermal quasiparticles n. The thermodynamical quanti- 
ties of interest can be deduced from the Helmholtz free 
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energy 



F{T) = -TlnZ(T), 



(1) 



where we express the temperature T in units of MeV and 
ks is set to unity. From this bridge equation connecting 
statistical mechanics and thermodynamics, we may cal- 
culate the entropy S, the average excitation energy (E), 
the heat capacity Cy, and the chemical potential /i at 
the temperature T by 



S{T) 
{E{T)) 
Cv{T) 



dF 

F + TS 
fd{E)\ 
V dT )y 
dF 
dn ' 



(2) 
(3) 

(4) 

(5) 



where n is the number of thermal quasiparticles. 

In the evaluation of the partition function we as- 
sume spin 1/2 fermions scattered into a doubly degener- 
ated single particle level scheme with equal energy spac- 
ing e. The level scheme has infinitely many levels, and 
it is one scheme for protons and one for neutrons, as in- 
dicated in Fig. |l|. These conditions simulate the Nilsson 
single particle scheme for deformed nuclei. 

We will first construct the partition functions z\ for 
n identical spin-up fermions placed in one and the same 
single particle scheme. We assume that n quasiparticles 
are available, with the core of Cooper pairs as a reser- 
voir. The total partition functions can then be expressed 
by these basic partition functions. The scattering of one 
spin-up fermion into the level scheme gives the partition 
function 



^exp(-*e/T) 



1 



1=0 



l-exp(-e/T)- 



(6) 



When two identical spin-up fermions are placed into the 
same level scheme, the summing is restricted due to the 
Pauli principle by 



zl = ^ exp(-ze/T) ^ exp(-j6/r) 

1=0 j=i+l 

^ 1 exp(-e/r) 

1 - exp(-2e/r) 1 - exp(-e/r) 

_ J exp(-e/r) 
U-exp(-2e/r)' 

Generally, for n identical spin-up fermions we find 

t exp(-(n-l)e/r) 



1 — exp(— ne/T) 



(7) 



(8) 



We now allow spin up and down for the fermions, and 
evaluate the corresponding partition function z„. For 



one fermion with spin up or down we simply obtain a 

T I 

degeneration of two, since z[ = z'{: 



Zl 



z{ = 2z\. 



(9) 



If two fermions occupy the same level scheme, the parti- 
tion function becomes more complicated. With the spins 
of the two fermions antiparallel (m = m\ -I- mi = 0), no 
Pauli blocking is present and we obtain 



zi(m = 0) = zjzj = 



(10) 



The contribution with parallel spin invoke the Pauli prin- 
ciple giving 



Ziijn = ±1) — z\ + Z2 



2z. 



ill) 



It is now straightforward to evaluate 2„ for n fermions, 
allowing spin up and down. For convenience we list the 
seven lowest partition functions here: 



Zl 

Z2 

zz 

Zi 
Zb 

z& 

Z7 



2z{ 
2zl 



(4)' 

2z\z\ 
2zr^ -f- 2z]^z^ 



2zl 
2zl 



2zl 



2zlz\ 



(4)' 

2Z3 Z2 
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{4)- 



(12) 



The partition function for nuclei with protons and neu- 
trons may now be constructed from the partition func- 
tions Zn- In Fig. the relevant configurations for the 
lowest quasiparticle excitations in the proton-odd system 
are illustrated. For this one-quasiparticle system, there 
is only one proton partition function to consider. For 
the three-quasiparticle case the ttu'^ and the tt^ partition 
functions have to be included. For the two-quasiparticle 
case of the even-even system we assume that either a tt or 
p pair can be broken. In the four-quasiparticle case there 
are n'^i/^, tt^ and z^"* configurations to take into account. 
From similar arguments for higher number of quasipar- 
ticles, we construct the seven lowest partition functions 
for proton-odd and even-even systems by 

Zl = Zl 

Zl = 2z2 

Zz = zz^ Z1Z2 

Zi = 2zi + {zif 

■^5 = ^5 + Z3Z2 + ZlZi 

Zq — 2zq + 2Z4Z2 

Z7 = Z7 + Z5Z2 + Z3Z4 + zizg. (13) 
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The partition functions for the neutron-odd system are 
identical to the functions for the proton-odd system. 

The partition functions for the odd-odd system can be 
constructed in a similar manner. The two-quasiparticle 
partition function is simply (zi)^, giving an entropy of 
twice the single quasiparticle entropy of the odd-mass 
system. Higher numbers of quasiparticles yield the fol- 
lowing partition functions for the odd-odd case: 

2z3Zi 
= 2z5Zi + (za)^ 
Zg = 2z7Zi + 2Z5Z3. (14) 

The single quasiparticle entropy is easily computed 
from our model with Eqs. (l,2Jl3) as 



5i(T) 



6/r 



In- 



exp(e/T) - 1 1 - exp(-e/T) 



(15) 



At higher excitation energies more quasiparticles have 
to be considered and the mathematical expressions be- 
come more complicated. In Fig. ^ the calculated n- 
quasiparticle entropy Sn is shown as function of e/T. 
The entropy is seen to scale rather well with the num- 
ber of quasiparticles, and there is a weak Pauli-blocking 
effect. The reason for this appearent extensivity is the 
onset of new degrees of freedom with increasing number 
of quasiparticles. In our formalism, this is manifested in 
the additional terms in the partition function appearing 
from the inclusion of spin up/down quasiparticles and 
the breaking of proton/neutron pairs. 

This approximate extensivity of the entropy can be uti- 
lized to estimate the many quasiparticle partition func- 
tion as a simple product of n quasiparticle partition func- 
tions Zn ~ {ZiY'-. The quasiparticle entropy is then 
roughly given by 



(16) 



where we argue in Sect. IV that the Si value is experi- 
mentally available for most nuclei. 

To create quasiparticles costs energy. In order to break 
one Cooper pair of the core the energy 2A is necessary. 
Taking this energy into account we can finally write the 
partition functions of even-even (ee) , odd (odd) and odd- 
odd (00) nuclei as weighted sums of the n-quasiparticle 
partition functions: 



(17) 





= 1 + 






^odd 


= Zi 






z°° 


= Z2 




+ e-^^/^Z, 



provided an ideal core of Cooper pairs, not changing with 
the number of quasiparticles. If the system is extensive, 
we may write the partition functions in closed form. By 
using the property of extensivity, Zn = (-^i)", and per- 
forming the sum up to N broken Cooper pairs^, we obtain 

l-(e-^/^Zi)^~+2 



1 - (e-^/^Zi)2 



{ZiY z- 



(18) 



The latter partition functions are good approximations 
for the expressions of Eq. (|l^) for up to ^2 broken 
pairs, however, for more than ~4 quasiparticles, the Pauli 
blocking cannot be neglected. 

Figure ^ shows the entropy as function of tempera- 
ture calculated from our model, allowing up to 5 broken 
nucleon pairs. Also calculations assuming extensivity of 
the entropy, in the sense Z„ ~ (^1)", are shown. There 
are serious discrepances between the non-extensive and 
extensive approaches for T > 0.3 MeV. 

There are limitations of the present model. The ex- 
cited states are assumed built by a few quasiparticles, 
while the other particles are forming Cooper pairs in the 
underlying core. This picture will break down for («) 
light nuclei, (ii) around closed shells and {Hi) at high 
excitation energy. The first two points rest on the fact 
that the single particle energy spacing should be small or 
comparable to the pairing interaction in order to build 
a superfluid phase. The last point is connected to the 
thermal quenching of pair correlations, where many va- 
lence particles participate in the excitation. A fingerprint 
of the breakdown of the thermal quasiparticle picture is 
when the entropy in the odd-mass and even-even systems 
becomes essentially^ equal. These limitations have been 
addressed in Ref. 



III. DETERMINATION OF LEVEL DENSITY 
ANCHOR POINTS 



The nuclear level density p is scarcely known as func- 
tion of excitation energy. Recently Q|, there has been 
published compilations of level density parameters. In 
our work we show that it is possible to determine two 
level density points quite reliably for numerous nuclei. 
To demonstrate the technique. Fig. ^ shows the two ex- 
tracted anchor points (filled data points) for ^''^Yb, to- 
gether with the level density deduced from known dis- 
crete levels (solid lines). The first point {Ei,pi) is based 



^For N ^ 00 broken Cooper pairs our simple model breaks 
down since the partition functions become infinite at the tem- 
perature e/T = - ln[l - 2 exp(-A/r)]. 
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on counting known discrete levels at low excitation en- 
ergy, giving an anchor point in the excitation energy re- 
gion of < < 2 MeV. Information on known levels is 
taken from the database of Ref . j|] , where we have ana- 
lyzed levels of ~ 1600 nuclei. The second point (i?2,P2) 
is estimated from the average neutron resonance spac- 
ing at the neutron binding energy _B„, where we use the 
compilation of Iljinov et al. 

Figure || also includes the data points ^ (open dr- 
ies) obtained with the method of Ref. 0. The dashed 
straight line describes the data rather well. This line will 
be referred to as the constant-temperature level-density 
formula, defined by 



piE) = Ce^/^ 



(19) 



with T-i = (lnp2 - lnpi)/(£'2 - Ei) and C ^ 
pi exp(— i?i/T). This line connecting the two anchor 
points is used to determine the nuclear temperature T. 

In the extraction procedure of {Ei, pi), an initial slope 
is predicted for the straight line; in the case of ^^^Yb, the 
inverse of T = 0.57 MeV. Then the two highest values of 
the count-based In p relative to this line is found^ . With 
maximums at the energy bin numbers a and 6, we find 
the lower anchor point by 



Pi = cxp 



1 



6-, 



(20) 



with energy Ei in bin number (a + b)/2. In the evalu- 
ation of the uncertainty of pi, the number of levels N 
(cTjv = V^) and the uncertainty of T (cft = O.IT) are 
taken into account. Nuclei far from /3-stability may have 
been poorly investigated due to low reaction cross sec- 
tions, and the extracted level density and uncertainties 
of these nuclei should be adopted with caution. 

The second point (i?2,P2) is estimated from the aver- 
age neutron resonance spacing at -B„. Assuming neutron 
£ — capture on a nucleus with target spin and parity 
levels with the same parity and two spins /± 1/2 are 
populated. In order to transfer the spacing D between 
these levels to the level density for all spins and parities, 
we use the spin distribution of Gilbert and Cameron B 



n -1 



le 2a 



(21) 



The spin cut-off parameter a is defined through cr^ = 
0.0888^2/3^^^ ^j^gj.g i-jjg back-shifted energy \s U = 

E— Ebs . The level density parameter a and the back-shift 
parameter E^s are defined by a = 0.21^° ®^ MeV^^ and 
Ebs = + A, respectively, where the back-shift correc- 
tion is given by Cbs = —6.6A~^-^'^ according to Ref. 



The pairing gap parameter is estimated by A = 12^ 
MeV. 

In Figs. we show the anchor points extracted for 
nine isotopes of germanium, gadolinium and uranium, re- 
spectively. The above procedure has been tested carefully 
and found reliable, as demonstrated in these figures. In 
Table 1 the two anchor points are listed together with the 
temperature T for the nuclei where both anchor points 
are known. We also see from Table 1 and Figs, that 
the extracted value of T is approximately equal for neigh- 
boring nuclei, thus, giving confidence to the methods ap- 
plied for establishing the two anchor points. 



IV. EXTRACTION OF ENTROPY 

The measured level density p{E) is proportional to the 
number of states accessible to the nuclear system at exci- 
tation energy E. Thus, the entropy in the microcanonical 
ensemble is given by 



S{E)^HpiE)/po)^liipiE) + So. 



(22) 



Applied to Eq. (p^), parameter T is interpreted as the 
thermodynamical temperature in the microcanonical en- 
semble, since by definition: 



{T{E)) 



{dS/dE)v 



(23) 



The constant shift parameter 5*0 = — In po can be ad- 
justed to fulfill the third law of thermodynamics; 5" ^ 
when r 0. The condition T = probably holds for 
the rotational ground band of deformed even-even nu- 
clei. However, also other interpretations for T = are 
possible, and we will not fix the Sq parameter to a spe- 
cific value in this work, but instead discuss the relative 
entropy S' — 5*0 = In p. 

The level density in between the two anchor points is 
not very well known. Recently [|l 10|, the level densities 
for i6i,i62j)y^ leegj. ^^^^ i7i,i72Yb have been measured. 

A common feature for these nuclei, is that the level den- 
sities develop more like the constant temperature level 
density formula of Eq. (|l^), than the back-shifted Fermi- 
gas level densit y as pointed out in the discussion of 
Fig. ^ and Ref. j| Ll|. This constant temperature behavior 
is interpreted |l|]10[| as the breaking of Cooper pairs with 
increasing excitation energy. 

In the following, we will discuss the microcanonical en- 
tropy for given excitation energy E. In order to establish 
a common play ground, we extrapolate the first and sec- 
ond level density anchor points to excitation energies 1 
and 7 MeV, using the constant temperature formula of 



^We pick the highest values since it is more probable that 
levels have escaped detection than the opposite. 



4 



Eq. (p^). Figures ^ and ^ display the entropy S — Sq 
based on the two anchor points evaluated at 1 and 7 
MeV, respectively. The data are plotted as function of 
the mass number, where nuclei along the /^-stability line 
are chosen. 



V. THERMAL QUASIPARTICLE PROPERTIES 
A. Regions of quasiparticle extensivity 

In the lower panel of figure ^ the difference in en- 
tropy between odd-mass and even-even systems AS* — 
S'odd — 'S'oo is shown. The data show fluctuations which 
are partly due to shell closures, incomplete knowledge 
of discrete levels at low excitation energy, and struc- 
tural changes between various isotopes and isotones. For 
the mid-shell regions the odd-mass nuclei have about 
AS* ~ 1.7 higher entropy compared to the even-even 
system. This feature holds for odd-even as well as for 
even-odd nuclei, however, neutrons carries slightly more 
entropy than protons. In Fig. ^ the data are sensitive to 
the level density extracted from neutron resonance level 
spacing. The fluctuations are less pronounced at 7 MeV 
because the level density is based on many quasiparti- 
cle states, smearing out structural effects. Also here, as 
indicated on the figure by lines, an entropy difference 
of AS" ~ 1.7 is found. Here, we should remember that 
these findings are based on a constant temperature level 
density extrapolation. If a Fermi gas level density is ap- 
plied, e.g. the model of Egidy et al. shown in Fig. ^, the 
extracted microcanonical entropy difference AS" will be 
somewhat modified; it will increase at 1 MeV and de- 
crease at 7 MeV. However, as demonstrated in Ref. 
the experimental value of AS* is rather constant for exci- 
tation energies above 1.5 MeV. 

It is amazing that AS is roughly the same both at 1 
and 7 MeV of excitation energy. This means that even 
with 4-6 quasiparticles excited at 7 MeV the last 
thermal quasiparticle still carries AS* ~ 1.7. Thus, the 
experimental entropy scales with the number of quasipar- 
ticles. This property is also accounted for in our model'^. 
The regions of quasiparticle extensivity can be identified 
in Figs. ^ and ^ as entropy-gaps between the entropies of 
odd-mass and even-even nuclei. The extensivity regions 
are found for mid-shell nuclei with mass number A > 40. 

In odd-odd nuclei, several neutron energy spacing data 
are known ||]. As seen in Fig. |9| (see lines), a proton and 
a neutron give an entropy excess of around 2AS ~ 3.4, 
as expected from the discussion in Sect. II. A similar 



study from known discrete levels is difficult, due to ex- 
perimental limitations in building level schemes with sev- 
eral hundreds of levels per MeV. However, in the case 
of the well studied odd-odd 1^2,154^^ isotopes the lower 
anchor point seems to have been determined with ac- 
ceptable confidence. Figure ^ shows that an additional 
neutron outside the odd-proton core gives an entropy dif- 
ference as high as 2.5. The extracted temperatures are 
all lying within a narrow window of T = 0.57 — 0.60 MeV. 

The feature that AS" is approximately independent of 
excitation energy E is shown in Table 1. Here, constant 
temperature T for neighboring mid-shell nuclei means 
that the slope of \np is equal and, thus, A^* is approxi- 
mately independent of E. Fig. |Tl| confirms the systemat- 
ics, where neighboring nuclei exhibit the same {dS/dE) 
value. This indicates that the level density, and hence 
the entropy, has the same functional form as function 
of excitation energy. However, the feature disappears at 
shell closures. 



B. Constancy of quasiparticle entropy 

Figure |l^ shows the entropy difference for a neutron 
outside its even-even core for various rare earth elements. 
Although various mid-shell isotopes reveal quasiparticle 
entropy that deviates from AS 1.7 it is still surprising 
how constant this quantity appears to be. However, as 
shown in Fig. |lj, the quasiparticle around the = 82 
shell closure carries low entropy due to the limited num- 
ber of quantuum states (multiplicity) available for the 
system. 

Our model gives a qualitative interpretation of the con- 
stancy of the quasiparticle entropy. With the assump- 
tions given, the model shows that the entropy Sn is a 
function of e/T, only. In order to deduce experimentally 
the value of e/T, we recognize that the single particle 
level spacing e is connected to the level density parame- 
ter a by PI 



7r\ , 7r2 

y (5p +5n) ~ y5 



TT 

37' 



(24) 



where the single-particle level-density parameters for pro- 
tons and neutrons {gp and (?„) are assumed to be approx- 
imately equal. Thus, we obtain using Eq. (p3|) 



f 



3a [dE 



(25) 



In Fig. the quantity {dS/dE)/a is plotted as function 
of mass number A, using the a values of Ref. Q and the 
data of Fig. ^ For increasing A, the extracted values 



^We should keep in mind that the model is formulated in 
the canonical ensemble, where a certain temperature T cor- 
responds to a wide excitation energy range. 
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are found in the range of 0.10 — 0.08, which corresponds 
to e/T = 0.33 - 0.26 according to Eq. (||). This varia- 
tion gives a minor effect on the predicted single particle 
entropy 5*1; Fig. ^ shows that the variation in e/T corre- 
sponds to an increase of 5*1 from 2.7 to 3.0. 

The experimentally extracted single quasiparticle en- 
tropy AS" shows no clear systematic trend as function of 
A within the uncertainties. Essentially, the local varia- 
tions due to structural changes between neighboring nu- 
clei are more pronounced than the weak variation in Si 
given by the model. 

The quasiparticle entropy AS* could in principle de- 
pend on the odd-mass ground state spin, since an orbital 
with high spin has more coupling combinations to the un- 
derlying rotating core. For excitation energies around the 
neutron binding energy i3„, the ground state orbital par- 
ticipates no more than other orbitals around the Fermi 
level, so no resemblence of the ground state spin is ex- 
pected. However, around the lower anchor point [Ei, pi), 
the ground state orbital could play a role. In Fig. |l4| we 
see no significant effect for rare earth nuclei, and the same 
result applies to other mass regions. This indicates that 
many orbitals participate to build up the level density in 
the lower excitation region of odd-mass nuclei. 

To close this subsection, we point out that Fig. ^indi- 
cates an average quasiparticle entropy AS* of 1.5 and 
~ 1.9 for protons and neutrons, respectively. The va- 
lence protons and neutrons in this figure occupy orbitals 
in the Z = 50 - 82 and iV = 82 - 126 shells, containing 
32 and 44 orbitals, respectively. Thus, the single parti- 
cle level energy spacing is expected to scale roughly as 
ep/e„ = 44/32 - 1.4. According to Fig. |, a 40% varia- 
tion in e/T around the value e/T — 0.3 gives a change in 
Si of about 0.4. Therefore, this effect could be the reason 
why the protons carry slightly less entropy compared to 
the neutrons. 



VI. APPLICATION OF THERMAL 
QUASIPARTICLES 

There are several applications for the concepts pre- 
sented here. In analogy with the description of rotational 
nuclei [|l3|, the Helmholtz free energy F plays the same 
role as the energy in the rotating frame E' . To visualize 
the analogy, the two cases read 

F = E-TS (26) 
E'^E-Loh, (27) 

where u is the rotational frequency and Ix is the spin 
alignment along the rotational x-axis. Generally, 5" and 
Ix are positive quantities, so that F and E' decreases 



*The core described by Zcoie is not necessary a cold core, 
but may already have pairs broken. 



with increasing T and w, respectively. A decoupled 
Cooper pair with high Ix may, at a certain critical fre- 
quency Wc, become more favoured than its coupled coun- 
terpart, giving rise to the rotational backbending phe- 
nomenon. 

In a similar manner, we will here demonstrate that the 
free energy of a system with two additional excited quasi- 
particles becomes favoured at a critical temperature Tc- 
This temperature depends on the pairing gap A and the 
quasiparticle entropy AS'. If both these quantities would 
be independent of the number of quasiparticles excited, 
a break-up of all Cooper pairs should happen at one and 
the same critical temperature Tc, giving rise to a strong 
quenching of the pair correlations. 

Let us assume that the partition function for an even- 
even nucleus is described by Z{T) = Zcotc{T)- This par- 
tition function represents the lowest free energy F{T) 
at low temperatures T. We now introduce the break- 
ing of one pair relative to this core^. Assuming the two 
quasiparticles to be independent of the core, we may de- 
scribe the excited system by a factorization of the parti- 
tion function into 

Z* w Zcoro^2e"'^/^, (28) 

where Z2 is the two quasiparticle partition function. 
With this assumption, the difference in Helmholtz free 
energy becomes 

Ai^ = -T{\tlZ* - InZcorc) = -rinZa + 2 A, (29) 

and factorizing Z2 assuming extensivity, gives 

AF = -TlnZi2 + 2A = 2(Fi-|-A), (30) 

where Fi is the free energy of one quasiparticle. This cor- 
responds to the chemical potential defined (see Eq. (||)) 
as the energy required to bring one quasiparticle out of 
the reservoir: 

^(T) AF/2 = Fi A. (31) 

We have applied this formalism to the experimental 
data of I6i,i62j3y^ j^gf_ r^^ie partition functions 
for each nucleus are calculated by 

Z{T) = J2 AEpiE)e-^/^ (32) 

E 

where AE is the experimental energy bin. The assump- 
tions made for evaluating Z is given elsewhere The 
corresponding Helmholtz free energy F for ^^^Dy is dis- 
played in the upper panel of Fig. |l^. The upper energy 
curve, based on data for ^^^Dy, is given by F(^^^Dy) + 
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A„ , where the pairing gap parameter is determined from 
neutron separation energies Sn by 

A„ = ^\SniN + 1, - 2Sn{N, Z) + 5„(7V - 1, Z)|, 

(33) 

giving A„ = 0.917 MeV for ^'^^Dy. Thus, the energy 
curve of ^^-'^Dy can be interpreted as the free energy for 
an even-even system with one extra quasiparticle. We 
therefore identify Fi + A„ in Eq. (pO) with the difference 
Fe6iDy) + A„-F(i62Dy). 

In the lower panel of Fig. |lj the chemical potential 
/i = i^i + A„ is shown as function of temperature. The 
higher entropy in ^^^Dy compared to ^^^Dy decreases the 
free energy relative to the even-even core with increasing 
temperature. When /i(T) ~ 0, thermal quasiparticles 
can be formed without consuming free energy. Thus, 
from Fig. ^ we estimate the critical temperature^ for the 
quenching of pair correlations to be Tc = 0.42 MeV. The 
uncertainty of this number is mainly due to the extensiv- 
ity assumed in the evaluation of AF. Furthermore, the 
adoption of a chemical potential /i = A at T = 0, rests on 
the assumption that the Fermi level coincides with single 
particle energies in both nuclei. Also the extrapolation 
of the experimental level density to higher excitation en- 
ergies, see Ref. gives systematical errors. From these 
considerations we estimate a 15 % error in the Tc value. 

The deduced value of Tc = 0.42(6) is considered to be 
common to the two dysprosium isotopes analyzed. Re- 
cently |], values of Tc = 0.52(4) and 0.49(4) MeV were 
found for ^^^Dy and ^^^Dy, respectively, using the canon- 
ical heat capacity as thermometer for the phase transi- 
tion. 

Finally, we point out a curious connection between the 
experimental quasiparticle entropy AS and the critical 
temperature Tc deduced from a Fermi gas model with 
pairing. If we assume a constant AS and a constant dif- 
ferences A between the E functions for the two curves of 
Fig. KB, we find 



VII. CONCLUSIONS 

We have analyzed the level densities of nuclei based on 
counting discrete levels and from neutron scattering data. 
It was possible to obtain level density anchor points for 
about 280 nuclei. The level densities have been trans- 
formed into microcanonical entropy, revealing the phase 
space of the system at fixed excitation energies. 

The single quasiparticle entropy has been extracted 
from the odd-even entropy differences. The quasipar- 
ticles not coupled in Cooper pairs, exhibit an entropy of 
AS* ~ 1.7 which is relatively independent of the pres- 
ence of other quasiparticles. This apparent extensivity is 
valid for mid-shell nuclei and for all ground state spins. A 
simple model for thermal quasiparticles in the canonical 
ensemble also supports the property of extensivity. 

The single quasiparticle entropy AS" for mid-shell nu- 
clei shows no clear mass dependency. This is explained 
in our model by the fact that AS only depends on the 
ratio between the single particle spacing e and the tem- 
perature T, which is found to be almost independent of 
the mass number for mid-shell nuclei. 

For thermal systems, the Helmholtz free energy is anal- 
ogous to the energy in the internal frame of rotating nu- 
clei. This opens for the study of various quantities and 
critical phenomena as function of temperature. We there- 
fore believe that the concept of single quasiparticle en- 
tropy may provide a fruitful basis for the understanding 
of hot nuclei. 

The simple model presented here is only adequate for 
mid-shell nuclei. Around closed shells one expects ex- 
citing effects from the increasing single particle energy 
spacings. This will also influence the entropy difference 
between odd-mass and even-even nuclei. Therefore, a 
statistical description of the transition to closed shells 
would be of great interest. 

We wish to acknowledge the support from the Norwe- 
gian Research Council (NFR). 



T — ^— A 



(34) 



giving Tc = 0.59A for AS" = 1.7. This estimate co- 
incides almost exactly with the theoretical expression 
Tc = 0.57A given in Refs. ||l|,|l|. The theoretical fac- 
tor 0.57 is independent on the mass number, consistent 
with the property of the experimental single quasiparticle 
entropy. 
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Table 1: Level density anchor points and temperatures. 



Nuclide 


Si(MeV) 


pi(MeV-i) 


£;2(MeV) 


P2 (MeV-i) 


T(MeV) 




2.70 


6( 2) 


6.60 


0.27E+02(0.70E+01) 


2.53(63) 




1.62 


6( 3) 


6.96 


0.38E+02(0.12E+02) 


3.04(90) 




0.78 


3( 2) 


7.33 


0.18E+02(0.57E+01) 


3.4(14) 


26Mg 


4.92 


12( 7) 


11.09 


0.76E+02(0.24E+02) 


3.3(12) 


2«A1 


1.32 


5( 2) 


7.73 


0.77E+02(0.21E+02) 


2.29(43) 


^«Si 


2.22 


4( 3) 


8.47 


0.54E+02(0.17E+02) 


2.40(71) 


32 p 


3.36 


9( 2) 


7.94 


0.91E+02(0.37E+02) 


2.00(43) 


33s 


2.46 


5(2) 


8.64 


0.82E+02(0.21E+02) 


2.20(41) 


34s 


2.70 


4( 3) 


11.42 


0.16E+03(0.62E+02) 


2.34(51) 


35s 


3.06 


7( 2) 


6.99 


0.60E+02(0.22E+02) 


1.84(44) 


36(-;i 


3.24 


12( 2) 


8.58 


0.37E+03(0.12E+03) 


1.56(17) 


38(^1 


1.86 


10( 5) 


6.11 


0.22E+03(0.69E+02) 


1.39(26) 


«Ar 


0.36 


3( 2) 


6.10 


0.14E+03(0.37E+02) 


1.52(31) 




2.58 


21( 6) 


7.80 


0.12E+04(0.16E+03) 


1.30(10) 


«K 


3.90 


34( 6) 


10.10 


0.26E+04(0.26E+03) 


1.44( 7) 


42^ 


1.44 


37(12) 


7.53 


0.11E+04(0.28E+03) 


1.81(23) 




3.12 


11(3) 


8.36 


0.34E+03(0.82E+02) 


1.54(17) 


«Ca 


2.88 


29( 7) 


7.93 


0.56E+03(0.11E+03) 


1.71(19) 


*^Ca 


3.90 


17( 5) 


11.13 


0.21E+04(0.33E+03) 


1.50(10) 


«Ca 


1.68 


10( 4) 


7.42 


0.48E+03(0.98E+02) 


1.49(18) 


«Ca 


3.00 


9( 3) 


5.15 


0.57E+02(0.20E+02) 


1.14(33) 




0.24 


17( 9) 


8.76 


0.43E+04(0.33E+03) 


1.54(16) 


47Ti 


2.70 


54(23) 


8.88 


0.11E+04(0.19E+03) 


2.07(32) 




4.02 


28( 6) 


11.63 


0.32E+04(0.82E+03) 


1.60(11) 


49TJ 


2.58 


19( 7) 


8.14 


0.91E+03(0.17E+03) 


1.44(15) 


SO^j 


3.42 


12( 4) 


10.94 


0.16E+04(0.39E+03) 


1.54(14) 


5lTi 


1.92 


6( 2) 


6.37 


0.24E+03(0.14E+03) 


1.23(23) 


Sly 


3.54 


59(23) 


11.05 


0.86E+04(0.20E+04) 


1.51(14) 


62y 


1.32 


9( 3) 


7.31 


0.13E+04(0.16E+03) 


1.22( 8) 


51 Cr 


1.14 


6( 2) 


9.26 


0.13E+04(0.22E+03) 


1.49(12) 


'«Cr 


1.86 


7( 2) 


7.94 


0.67E+03(0.13E+03) 


1.32(10) 


s^Cr 


2.64 


8( 3) 


9.72 


0.11E+04(0.38E+03) 


1.46(15) 


55 Cr 


0.42 


6( 3) 


6.25 


0.30E+03(0.54E+02) 


1.47(23) 


5«Mn 


0.66 


15( 5) 


7.27 


0.20E+04(0.32E+03) 


1.35(11) 


sspe 


0.66 


3( 2) 


9.30 


0.11E+04(0.22E+03) 


1.48(19) 




0.66 


4( 2) 


7.65 


0.83E+03(0.15E+03) 


1.33(11) 




2.46 


10( 4) 


10.05 


0.19E+04(0.52E+03) 


1.45(13) 




0.84 


8(4) 


6.58 


0.57E+03(0.25E+03) 


1.36(20) 




1.08 


19( 3) 


7.49 


0.42E+04(0.61E+03) 


1.19( 5) 


«iCo 


1.44 


12( 4) 


9.32 


0.22E+05(0.16E+04) 


1.06( 5) 




1.02 


6( 2) 


9.00 


0.17E+04(0.29E+03) 


1.40( 8) 




3.84 


30( 8) 


11.40 


0.51E+04(0.12E+04) 


1.47(10) 




2.04 


32(11) 


7.82 


0.14E+04(0.26E+03) 


1.53(16) 


62 Ni 


3.66 


27( 6) 


10.60 


0.51E+04(0.84E+03) 


1.32( 7) 


63 Ni 


1.14 


9( 5) 


6.84 


0.11E+04(0.24E+03) 


1.18(13) 




0.48 


4( 3) 


6.10 


0.75E+03(0.17E+03) 


1.05(15) 


«^Cu 


1.68 


49 ( 7) 


7.92 


0.11E+05(0.14E+04) 


1.16( 4) 


««Cu 


0.30 


13( 6) 


7.07 


0.74E+04(0.16E+04) 


1.07( 8) 


s^Zn 


0.96 


25(10) 


7.98 


0.67E+04(0.20E+04) 


1.26(11) 


«^Zn 


1.86 


31( 7) 


7.05 


0.36E+04(0.66E+03) 


1.09( 6) 


««Zn 


3.42 


37( 8) 


10.20 


0.21E+05(0.26E+04) 


1.07( 4) 


«9Zn 


1.20 


18( 6) 


6.48 


0.27E+04(0.46E+03) 


1.06( 8) 


^iZn 


0.60 


8( 3) 


5.83 


0.33E+04(0.58E+03) 


0.86( 7) 


™Ga 


1.14 


57(12) 


7.65 


0.35E+05(0.88E+04) 


1.02( 5) 


^^Ga 


0.54 


44( 7) 


6.52 


0.21E+05(0.43E+04) 


0.97( 4) 


^iGe 


1.32 


51(14) 


7.42 


0.13E+05(0.53E+04) 


1.10(10) 


^^Ge 


0.84 


42(13) 


6.78 


0.12E+05(0.20E+04) 


1.06( 7) 
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2.64 


28( 5) 


10.20 


0.74E+05(0.80E+04) 


0.96( 3) 




0.24 


20( 9) 


6.51 


0.66E+04(0.14E+04) 


1.09( 9) 


■^^Gc 


0.60 


20( 9) 


6.07 


0.32E+04(0.44E+03) 


1.08(10) 


^®As 


0.54 


82(15) 


7.33 


0.89E+05(0.12E+05) 


0.97( 3) 




1.08 


62(16) 


8.03 


0.67E+05(0.39E+05) 


1.00( 9) 


"Se 


0.96 


34( 8) 


7.42 


0.36E+05(0.80E+04) 


0.93( 4) 




2.52 


29 ( 7) 


10.50 


0.15E+06(0.30E+05) 


0.94( 4) 


™Se 


0.66 


24( 7) 


6.96 


0.16E+05(0.33E+04) 


0.97( 6) 


®^Se 


0.60 


12( 7) 


6.70 


0.81E+04(0.35E+04) 


0.93(10) 


*^Se 


0.48 


17( 8) 


5.82 


0.39E+04(0.73E+03) 


0.98(10) 


*°Br 


0.36 


75(18) 


7.89 


0.17E+06(0.21E+05) 


0.98( 3) 


^^Br 


0.90 


43(11) 


7.59 


0.98E+05(0.18E+05) 


0.87( 4) 


™Kr 


0.84 


44( 9) 


8.36 


0.13E+06(0.37E+05) 


0.94( 4) 


SI T ^ 


0.96 


41(14) 


7.88 


0.11E+06(0.12E+05) 


0.88( 4) 


^^Kr 


0.84 


11( 3) 


7.46 


0.80E+05(0.51E+05) 


0.75( 6) 


*®Rb 


1.08 


32(11) 


8.65 


0.77E+05(0.16E+05) 


0.97( 5) 


®*Rb 


0.30 


12( 7) 


6.08 


0.57E+04(0.16E+04) 


0.94(10) 


®^Sr 


1.14 


16( 4) 


8.53 


0.83E+05(0.32E+05) 


0.86( 5) 




1.50 


9( 4) 


8.43 


0.34E+05(0.11E+05) 


0.85( 6) 




2.52 


16( 9) 


11.11 


0.60E+05(0.67E+04) 


1.04( 7) 




2.04 


20( 9) 


6.37 


0.25E+04(0.49E+03) 


0.90( 9) 


90 Y 


0.96 


9( 4) 


6.86 


0.15E+05(0.28E+04) 


0.80( 5) 




2.28 


37(15) 


7.20 


0.52E+04(0.10E+04) 


1.00( 9) 




1.44 


8( 6) 


8.64 


0.24E+05(0.68E+04) 


0.90( 9) 


®^Zr 


1.80 


16(4) 


6.73 


0.13E+05(0.36E+04) 


0.74( 4) 


^^Zr 


1.80 


18( 6) 


6.46 


0.13E+05(0.25E+04) 


0.71( 4) 


''''Zr 


1.20 


8( 6) 


5.58 


0.11E+05(0.30E+04) 


0.61( 7) 


®*Nb 


1.14 


61(13) 


7.23 


0.98E+05(0.14E+05) 


0.82( 3) 


^^Mo 


2.76 


103(21) 


8.07 


0.20E+05(0.38E+04) 


1.01( 5) 


^5 Mo 


0.90 


17( 9) 


7.37 


0.30E+05(0.97E+04) 


0.87( 7) 


''•'Mo 


2.58 


44(11) 


9.15 


0.86E+05(0.12E+05) 


0.87( 3) 


''''Mo 


0.90 


20( 5) 


6.82 


0.35E+05(0.61E+04) 


0.80( 3) 


^*Mo 


2.46 


51(12) 


8.64 


0.99E+05(0.15E+05) 


0.82( 3) 


''"Mo 


0.78 


35(10) 


5.93 


0.34E+05(0.84E+04) 


0.75( 4) 




0.42 


44(11) 


5.40 


0.46E+05(0.56E+04) 


0.72( 3) 


100 


0.42 


75(14) 


6.76 


0.28E+06(0.55E+05) 


0.77( 3) 


^°°Ru 


2.34 


46( 9) 


9.67 


0.24E+06(0.53E+05) 


0.86( 3) 


^°^Ru 


2.16 


24( 4) 


9.22 


0.46E+06(0.61E+05) 


0.72( 2) 


^°^Ru 


0.36 


34( 8) 


6.23 


0.64E+05(0.19E+05) 


0.78( 4) 


^°*Ru 


1.38 


11(6) 


8.91 


0.11E+07(0.56E+06) 


0.65( 4) 


"'■'Rh 


0.48 


121(22) 


7.00 


0.56E+06(0.81E+05) 


0.77( 2) 


105 


0.66 


33( 6) 


7.09 


0.16E+06(0.25E+05) 


0.76( 2) 


106 


2.52 


56(11) 


9.56 


0.65E+06(0.97E+05) 


0.75( 2) 




0.54 


24(11) 


6.54 


0.14E+06(0.48E+05) 


0.69( 5) 


108 


2.16 


23( 6) 


9.22 


0.72E+06(0.13E+06) 


0.68( 2) 


109 p^ 


0.48 


46(14) 


6.15 


0.19E+06(0.77E+05) 


0.68( 4) 




0.36 


28 ( 7) 


5.75 


0.10E+06(0.17E+05) 


0.66( 3) 


108 


0.60 


115(22) 


7.27 


0.70E+06(0.15E+06) 


0.77( 3) 


""Ag 


0.48 


87(19) 


6.81 


0.98E+06(0.20E+06) 


0.68( 2) 


^°''Cd 


0.66 


14( 4) 


7.93 


0.31E+06(0.85E+05) 


0.73( 3) 


^°"Cd 


0.96 


26 ( 7) 


7.35 


0.34E+06(0.91E+05) 


0.67( 3) 


"^Cd 


0.84 


38( 9) 


6.98 


0.26E+06(0.42E+05) 


0.70( 2) 


"^Cd 


2.88 


190(57) 


9.40 


0.86E+06(0.16E+06) 


0.77( 3) 


"^Cd 


1.08 


87(19) 


6.54 


0.20E+06(0.44E+05) 


0.71( 3) 


"*Cd 


1.86 


19(10) 


9.04 


0.83E+06(0.13E+06) 


0.67( 3) 


Cd 


0.60 


20( 6) 


6.14 


0.17E+06(0.30E+05) 


0.61( 2) 


"^Cd 


0.24 


24(10) 


5.77 


0.10E+06(0.25E+05) 


0.66( 4) 




1.08 


68(16) 


7.28 


0.10E+07(0.15E+06) 


0.65( 2) 


"«In 


0.90 


80(15) 


6.78 


0.73E+06(0.52E+05) 


0.65( 1) 




1.14 


19( 7) 


7.75 


0.28E+06(0.95E+05) 


0.69( 4) 
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0.60 


12( 7) 


7.55 


0.19E+06(0.68E+05) 


0.72( 5) 


"^Sn 


2.22 


16( 7) 


9.56 


0.54E+06(0.27E+06) 


0.70( 5) 


"^Sn 


1.56 


44(17) 


6.95 


0.11E+06(0.40E+05) 


0.69( 5) 


"®Sn 


2.22 


27(11) 


9.33 


0.53E+06(0.13E+06) 


0.72( 3) 


"^Sn 


1.08 


18( 6) 


6.49 


0.70E+05(0.19E+05) 


0.66( 3) 




1.74 


12(4) 


9.11 


0.38E+06(0.83E+05) 


0.71( 3) 




1.02 


24( 8) 


6.17 


0.42E+05(0.13E+05) 


0.69( 4) 




1.02 


28(13) 


5.95 


0.26E+05(0.13E+05) 


0.72( 7) 




0.96 


20( 9) 


5.73 


0.17E+05(0.82E+04) 


0.71( 7) 




0.24 


62(16) 


6.81 


0.53E+06(0.72E+05) 


0.73( 2) 




0.48 


61(16) 


6.47 


0.34E+06(0.41E+05) 


0.69( 2) 


123 rpg 


0.60 


25(10) 


6.93 


0.32E+06(0.12E+06) 


0.67( 4) 


124 rpg 


2.76 


112(22) 


9.43 


0.14E+07(0.34E+06) 


0.71( 2) 


125 rpg 


0.54 


24( 8) 


6.57 


0.20E+06(0.50E+05) 


0.67( 3) 


126 rpg 


2.28 


46(10) 


9.11 


0.66E+06(0.83E+05) 


0.71( 2) 


127 rpg 


0.72 


29(13) 


6.29 


0.74E+05(0.14E+05) 


0.71( 4) 


129rpg 


0.60 


16( 5) 


6.09 


0.22E+05(0.70E+04) 


0.75( 5) 


131 rpg 


0.96 


25(10) 


5.93 


0.79E+04(0.13E+04) 


0.86( 7) 


128j 


0.30 


75(14) 


6.83 


0.52E+06(0.14E+06) 


0.74( 3) 


130 J 


0.24 


47(14) 


6.46 


0.32E+06(0.62E+05) 


0.71( 3) 




0.72 


45 ( 9) 


7.60 


0.76E+06(0.83E+05) 


0.71( 2) 




0.42 
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FIG. 1. Illustration of thermal excitations in an odd proton nucleus, where the horizontal bars indicate how much the orbitals 

arc filled on the average. In the case of no broken pairs (left diagram), the single proton has maximum number of quantuum 
states (multiplicity) available. The three quasiparticle excitation has to include either a Trf^ or a tC' configuration (two right 
diagrams) . Here, the Pauli blocking reduces the multiplicity and gives a lower quasiparticle entropy. 
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FIG. 2. The entropy 5„ of n quasiparticles. The entropy increases with increasing temperature T for fixed level spacing e. 
The lower part demonstrates that the Pauli blocking reduces the entropy per quasiparticle for n > 2. 
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FIG. 3. The total entropy calculated for even-even, odd-mass and odd-odd nuclei allowing a maximum of five broken pairs. 
The extensive and non-extensive calculations are performed with Eqs. (^sj) and ([l7|), respectively. The parameters used are A 
= 0.9 MeV and e = 0.15 MeV. 
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FIG. 4. Level densities in Yb, showing the two anchor points (filled circles). The solid line is based on the counting 
of discrete known levels. Also shown are the backshifted level density prediction of Egidy et al. [9] (dashed-dotted) and the 
constant temperature formula (dashed). 
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FIG. 5. The level density anchor points extracted for ®* '^^Ge. The lower point is based on known discrete levels. In cases 
where the second anchor point is known, the extracted temperature T is quoted. 
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FIG. 6. Anchor points for ^^^-woq^j^ ^^^^ pjg || 
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FIG. 7. Anchor points for 231-239^^ gee text of Fig. ||. 
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FIG. 8. Entropy as function of mass number at 1 MeV of excitation energy. The upper panel shows that the odd-even (filled 

triangles) and even-odd (filled squares) nuclei exhibit higher entropy than the cvon-cvon (open circles) nuclei. In the lower 
panel, the entropy difference between odd-mass and cvcn-ovcn nuclei is shown, giving an average of AS ^ 1.7 (solid line) for 
mid-shell nuclei. The appoximate locations of the neutron and proton magic numbers arc indicated. The data are based on the 
anchor points discussed in the text, and each data point is typically the average of 2-4 nuclei centered around the /3-stability 
line. 



21 



25 



22.5 



20 



17.5 



cn 15 
CO 

^ 12.5 

o 



10 



7.5 



2.5 



E 



eV 



N = 82 



Z = 50 



N = 50 



N = 28 
N = 20 



. ^ c 




i, . i 
i 

^0 °3 



Ml, 



.CD *o 







Ik A o 



N=126, Z = 82 



25 50 75 100 125 150 175 200 225 250 

Mass number A 



FIG. 9. Entropy as function of mass number at 7 MeV of excitation energy. The data are plotted for odd-odd (open triangles), 
odd-even (filled triangles), even-odd (filled squares) and even-even (open circles) nuclei, see text of Fig. ^. The three lines are 
drawn to indicate the entropy gaps between even-even, odd-mass and odd-odd nuclei. 
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FIG. 10. Level densities and entropy difference for odd-even 

i53.i55g^ and odd-odd is2,i54g^ nuclei. 
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FIG. 11. The derivative of entropy with respect to excitation energy as function of mass number. The derivatives are 
determined from the anchor points of nuclei along the line of /3-stability. The data are plotted for odd-even (filled triangles), 
even-odd (filled squares) and even-even (open circles) nuclei. 
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FIG. 12. The entropies extracted at 1 MeV of excitation energy for odd-even (filled triangles) and even-even (open dries) 
rare earth isotopes. The data demonstrates that the entropy difference A5 between odd-even and even-even isotopes is strongly 
reduced at the N = 82 shell closure. 
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FIG. 13. The experimentally available quantity {dS/dE)/a is proportional to e/T, which determines the single quasiparticle 
entropy in our model. 
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FIG. 14. The entropy difference AS for odd-even (filled triangles) and even-odd (filled squares) nuclei relative to even-even 
nuclei as function of ground state spin. The data are taken for nuclei with Z = 60 — 78 and N = 90 — 110. The solid lines are 
the averages for protons (1.5) and neutrons (1.9). 



27 



1 000 
500 


^ -000 
c 

CD -1000 

CD 

^-1500 
-2000 



> 

CD 











\. F( 






F( 




~ 1 1 1 1 


1 1 1 1 





0.1 



0.2 



0.3 



0.4 



0.5 



> 

(U 

o 

CP 

o 

O 

o 

CD 



1000 



800 



600 



400 



200 



A = 0.92 MeV 



Tc = 0.42 MeV 



0.2 0.3 0.4 0.5 

"emperature T [MeV] 



FIG. 15. The experimentally deduced Helmholtz free energy for ^®^'^®^Dy extracted in the canonical ensemble. The critical 
temperature for the quenching of pair correlations is found at Tc = 0.42 MeV, where the chemical potential is /x ~ 0. 
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